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1. Introduction 

The Einstein-Maxwell system of equations plays a central role in relativistic astrophysics 
when describing spherically symmetric gravitational fields in static manifolds. In these 
situations we are modelling charged compact objects with strong gravitational fields such 
as dense neutron stars. Some recent investigations, including the models of Ivanov [1], 
Sharma et al |2j and Thirukkanesh and Maharaj [3], indicate that the electromagnetic 
field significantly affects physical quantities in relativistic stellar systems: equations 
of state, redshifts, luminosities, stability and maximum masses of compact relativistic 
stars. The presence of electric charge is a necessary ingredient in the structure and 
gravitational evolution of stars composed of quark matter as illustrated in the treatments 
of Mak and Harko [1] and Komathiraj and Maharaj [5] . These applications emphasise 
the need for exact solutions of the Einstein-Maxwell system in modelling the formation 
and evolution of charged astrophysical objects. Lasky and Lun |6l[7], amongst others, 
studied the role of electromagnetic fields in gravitational collapse, formation of black 
holes and the existence of naked singularities. Electric fields cannot be ignored in 
spherical gravitational collapse with phases of intense dynamical activity and particle 
interaction as shown in the treatments of Di Prico et al [8] and Herrera et al [9]. 
Maxwell's equations also play an important role in cosmological models in higher 
dimensions, brane world models and wormhole configurations. In this context see the 
treatment of De Felice and Ringeval [10] . 

When solving the Einstein field equations with neutral matter distributions, we 
often make the assumption that the spacetime is shear-free and spherically symmetric. 
Kustaanheimo and Qvist |Tl] were the first to present a general class of solutions. The 
generalisation to include the presence of the electromagnetic field may be easily achieved. 
The field equations are reducible to a single partial differential equation. A review of 
known charged solutions, with a Friedmann limit, is given by Krasinski [12]. Srivastava 
^13j and Sussmann [HI [15] undertook a detailed study of the mathematical and physical 
features of the Einstein-Maxwell system in spherical symmetry. Wafo Soh and Mahomed 
[T6] used symmetry methods to systematically study the underlying partial differential 
equation. They showed that all previously known solutions can be related to a Noether 
point symmetry 

The main aim of this paper is to study the integrability properties of the underlying 
partial differential equation for the Einstein-Maxwell system using symmetry methods. 
Both Noether and Lie point symmetries of the governing equation are considered. 
Noether symmetries have the interesting property of being associated with physically 
relevant conservation laws in a direct manner via the well-known Noether theorem. 
The Lie symmetries are more general, providing a larger set of symmetry generators in 
general, but do not guarantee integrability and reduction to quadrature in a straight 
forward manner. In section 2, we reduce the Einstein-Maxwell field equations to a 
single nonlinear second order partial differential equation that governs the behaviour 
of charged fluids. This is achieved by utilising the generalised transformation due to 
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Faulkes |17j. The resulting partial differential equation can be treated as an ordinary 
differential equation as in the case of uncharged fluids. In section 3, we present a 
first integral of the governing equation obtained earlier by generalising the technique 
of Srivastava [13] first used for uncharged fiuids. This first integral is subject to two 
integrability conditions expressed as nonlinear integral equations. We transform the 
integral conditions into a fourth order differential equation. In section 4, we analyse the 
governing equation for Noether symmetries via its Lagrangian and this analysis yields 
a general Noether first integral for this equation. We then establish the relationship 
between the Noether first integral and the first integral (which had been obtained earlier 
using an ad hoc approach). In section 5, we undertake a comprehensive Lie symmetry 
analysis of the governing equation to investigate the conditions under which it can be 
reduced to quadratures. We show how the Noether results are a subset of the Lie 
analysis results. Lastly, in section 6, we discuss the results obtained and relate some 
invariant solutions to known results. 



2. Field equations 

We analyse the shear-free motion of a fiuid distribution in the presence of an electric 
field. It is possible to choose coordinates = (t, r, 9, 0) such that the line element can 
be written in the form 

which is simultaneously comoving and isotropic. Since an electromagnetic field is 
present, the Einstein field equations are supplemented with Maxwell's equations to 
describe a self-gravitating charged fiuid. The Einstein field equations for a charged 
perfect fiuid can be written as the system 

_ A2 1 / 4A,\ 
P - 3^ - ^ ( 2A^^ + -^r + — I - -^x^ 



— (-3A2 - 2Xu + 2UtXt) + ^ Urr + + - + — + A, 





(2 a) 




{2b) 




E^ 


(2c) 


^4g4A 




{2d) 




(3) 



p 

= UrXt - {Xt)r 

Maxwell's equations imply that 

E = r^e^-''*,., Er = ar^e^^ 

In the above p is the energy density, p is the isotropic pressure, a is the proper 
charge density and E is the total charge contained within the sphere of radius r 
centered around the origin. The matter variables are measured relative to the four- 
velocity fi"" = (e"'^, 0, 0, 0). Subscripts refer to partial derivatives. Observe that 
$j. = FiQ {= —Fqi) is the only nonzero component of the Faraday electromagnetic field 
tensor Fab = (t>b;a — (pa;b where we have chosen (pa = ($(t, r), 0, 0, 0) as the simplest 
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choice that allows for a nonvanishing electric field. The coupled equations f l2fflp -(l2( 
(|3]) comprise the Einstein-Maxwell system for the metric ([1]) in terms of the variables 
(j),p,E,a, I/, A. 

The Einstein-Maxwell field equations can be reduced to a simpler system as shown 
by Srivastava [18]. After considerable simplification the Einstein field equations (l2aj) - 
f l2dj) can be written in the equivalent form 

■"K + A? + i^)--|l, (4a) 



p = 




p = 


1 






= Xte- 


^1 






g3A+2/i{t) _ E 6 



^4 



J t 



2EV 

^4 



m 

(4c) 
(Ad) 



for a charged relativistic fiuid, where h = h{t) and F = F{r) are arbitrary constants of 
integration. We have the condition At 7^ from ( l46|) so that the model must be time 
dependent, and the spacetime cannot become static. Observe that fj4dj) is essentially 
the generalised condition of pressure isotropy which includes a term corresponding to 
the electric field. To find an exact solution to the system (14 fflp -( l4dl) we need to explicitly 
integrate fj4dj) to determine A, and cr follows from ([3]). 

To continue we need to express ( l4dj) in a simpler form. We use a transformation 
first utilised by Faulkes fT7] for neutral fiuids. We let 

X = (5 a) 

y = e-' {5b) 

g{x) = {5d) 

With the help of (15 fflP -f l5dj) we can write ( |4^ in a more compact form 

y" = f{x)y^ + g{x)y'^ (6) 

where primes denote differentiation with respect to the variable x. Equation ([6]) is the 
fundamental nonlinear partial differential equation which determines the behaviour of 
the self-gravitating charged fluid in general relativity. However, we can treat it as an 
ordinary differential equation as only derivatives with respect to x appear. U g = then 
we regain a neutral fluid which has been studied by Maharaj et al [19j and Wafo Soh 
and Mahomed [20], amongst others. 



3. A charged first integral 



It is possible to flnd a flrst integral of (|6]) without choosing explicit forms of the functions 
/(x) and g{x). We use integration by parts, an approach adopted by Maharaj et al [19] 
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in investigating the integrability properties of the field equation y" = f{x)y'^ of a neutral 
spherically symmetric shear-free fiuid. We can integrate ([6]) by parts to obtain 

y' = fiy' + 9iy' - Vnyy' + 2fniy'^ + 2 + 2 {gfjj), y' 

-K^y^-K,y^ + r^{t) (7) 

For convenience we have used the notation 

f{x)dx = fi / g{x)dx = gj (8) 



The result ([7]) is subject to the integrability conditions 

Ko = ^ffni + 2{ffu)j + gi {9a) 
Ki = gfin + 2{gfn), {9b) 

where Ko and Ki are constants and To{t) is an arbitrary function of integration. 
The integral equations fl9aj) - fl96j) may be simplified. On setting 

fin = a (10) 
we can rewrite (|96|) as 

ga + 2 {ga') j = Ki 
which has the solution 

g = 9oa-^ (11) 
where a = a{x) and go is an arbitrary constant. Equation flQoj) now becomes 

^a"'a + 2{a"'a')j+[goa-')^ = Ko 

or, if we want a purely differential equation, 

aa(™) + ^aV" = -^goa^'' (12) 

We have therefore established that the integrability conditions (I9fflj)-( l9^ can be 
transformed into the fourth order ordinary differential equation ( fT2|) . Using (l9a|l -( l96jl . 
( ITOl) and ( JTTl) we rewrite the first integral ([7]) of (E]) as 

4 

To{t) = y' - a"y^ + 2a yy' - 2ay'^ + -aa"y^ + goa'^y^, (13) 

subject to condition (fT2|) . A solution of (fT2|) will give a, and then / and g will be found 
from (fTOll and ( ITTl) respectively. A full analysis of this case can be found in Kweyama 
et al pT], and we will not repeat those results here. 
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4. Noether symmetries and integration 

Given that Noether's theorem generates first integrals in a direct manner, we now 
investigate (j6]) for Noether symmetries in order to generate first integrals. Note that, if 
a second order ordinary differential equation 

y" = N{x,y,y') (14) 
has a Lagrangian C{x,y,y') then ( IT^ is equivalent to the Euler-Lagrange equation 
d_ /dC\ _dC_^ 

dx \dy' J dy 
The determining equation for a Noether point symmetry 

X = ^{x,y)-^ + r]{x,y)-^ (15) 

corresponding to a Lagrangian £(x, y, y') of ffT^ is 

dx / dx 



where 

\dx dx J oy 

and F = F{x,y) is a gauge function. It is also known that if there is a Noether 
symmetry corresponding to a Lagrangian of an equation f|T^ . then f|T^ can be reduced 
to quadratures. This is a critical advantage of Noether point symmetries as emphasised 
by Wafo Soh and Mahomed [TB]. The Noether first integral 1^ associated with the 
Noether point symmetry (ITS!) is given by 

In = y)C + ir]ix, y) - y'^x, y)) Cy, - F (16) 

in terms of C and F. 

For equation ([6]) a Lagrangian is 

C = \y'' + \f{x)y' + \g{x)y' (17) 
The Lagrangian C admits the Noether point symmetry 



provided 



fe=^a' (19a) 

g = g,a~' {19b) 

f = a-'/' (/i - 3(71 J ca-'/'dx^ (19c) 

a'" = 4a-^/^c (^/i - 3^1 J ca'^^^dx^ (19 d) 

c = Co + Cix (19e) 

F=^a"y' + c'y (19/) 
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where a = a{x), b = b{x), c = c{x) and /i and gi are constants. From fll9cp and (119dp 
we have 

n'" 

f= — 
■' 4c 

When we differentiate fll9(jj) once with respect to x we obtain 

caa^'"^ + ^ca'a'" - ac'a'" = -Ug^c^a'^ (20) 

Applying the transformation 

X = A = ^'^i (21) 

Co + C,x' (Co + Cix)' ^ ^ 

reduces (120|) to 

^^(™) + ^A'A'" = -UgiClA-^ (22) 

which is in the form of f|T2l) . 

From flTB]) we have that (jH]) admits the Noether first integral 
■1 ,o 1 



In = a 



+ \y'a'^"'^ i^f, - 3g, J ca-(^/^)dx) + \g,a-Y 



2' 

+ ( la'y + c - ay') y' - -a"y^ - c'y (23) 



2 - " y " 4 

Substituting ffTOdl) in ([23]) yields 



In = cy - -a y +-ayy - -ay + j^y + -g^a y - cy. (24) 



On comparing our first integral (lT3ll and the Noether first integral ( 1241) we observe that, 
with c = 1/4, /o = /i and go = gu 

ro{t) = 41 N 

As a result, our earlier ad hoc approach yielded a first integral which is a special case of 
that obtained via Noether's theorem. This further implies that ( !T3|) admits the Noether 
symmetry 

F = 4a|- + {2a'y + 1) |- 
ox oy 

thus supporting the results of our approach, and the fact that ^ could be reduced to 
quadratures. 



5. Lie analysis 

As a final attempt at solving (jS]) we undertake a Lie symmetry analysis [22]- Lie 
symmetries are usually a larger set of symmetries for a problem as opposed to the set 
of Noether symmetries. However, the disadvantage is that no simple formula exists to 
find first integrals associated with Lie symmetries - direct integration of (often difficult) 
equations is usually needed. However, as we show here, this approach allows for a larger 
class of solutions than the Noether approach. 
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It is a simple matter to verify that (1181) is a symmetry of (jSD and the relationship 
among the functions a{x),b{x),c{x), f{x) and g{x) is given by the following system of 
ordinary differential equations 

a" = 2b' (25 a) 



b" = 2/c 



c" = 



af' + {2a' + b)f = -3cg 
ag' + {2a' + 2b)g = 



(255) 
(25 c) 
{25 d) 
(25 e) 



From fl25ai) 



2b = a' + a 

where a is an arbitrary constant, from (125 a\i and fl25fe|) 

a'" 
f = — 

^ Ac' 

from (125 c\i we have (I19ep again, and finally from (125 ep 

adx\ 



g = g2a exp 



(26) 



(27) 



where g2 is an arbitrary constant. By integrating ( 125 J) and using (1271) we obtain an 
alternative form of fl26|) . and it is 



/ 



a exp I - 



It) /.-3./..-<-exp(/^).. 



(28) 



where /2 is an arbitrary constant of integration. From fl25dj) . fl26|) and fl27|) we observe 
that a is a solution of the equation 



caa^'""^ + 



' ba' 



da 



^I2g2c'a-^exp(^- (29) 



Alternatively (1291) may be obtained by equating equations (!26l) and (!28|) and 
differentiating once with respect to x. 

Observe the similarity between the results obtained here in the case of Lie analysis 
and those obtained in the previous section on Noether symmetries. The main difference 
is the occurrence of the parameter a in the Lie analysis. Thus, we have a more general 
result (as expected) in the Lie analysis as compared to the Noether analysis. In addition, 
this result is even more general than that obtained in Kweyama et al [21] as our c(x) is 
nonconstant. 

d 

The transformation which converts the symmetry ( ITSl) to —— makes ^ 
autonomous. The simplest expression of this transformation is 



X 



dx 
a 



Y 



y exp - 



bdx^ 



exp 



bdx 



dx 



(30) 



Using f l5U]) equation (E]) becomes 



a 



Y" + aY'+{M+ — ]Y = f^Y^ + g^Y^ + 



(31) 
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where /2 and g2 are arbitrary constants introduced in (P7|) and (1251) respectively. The 
quantities M and are arbitrary constants that arise in integrations of f l29|) and are 
given by 

M = -aa" - ]a'^ - 2/2/ + 3g2p (32) 



and 



2 4 



-1/2 / / Cidx\ ( I \ , 1 

A/ = —a ' exp — / \ac a c H — ac 

^ ' ' 2a \ 2 2 , 



where 



iaa" - + x) ^ + ~ ^^^^^ ^^^^ 

ca-'^^/'^ exp J (34) 

Note that in the neutral perfect fluid case we must have g2 = N = and fl^ reduces 
to that of Maharaj et al [19] . However, it is difficult to make direct comparison to the 
results therein as the full equations are not always given, and the equation referencing 
is not always clear. 

To proceed further, we need to analyse equation fl29|) . We note that the form of 
( 129|) was obtained under the assumption that c 7^ 0. We consider both nonzero c and 
vanishing c in turn in our subsequent analysis. 

5.1. Case I: c ^ 

When Co 7^ and Ci = 0, (l29l) may be written as 

aa(™) + i (5a' + a) a'" = -12g2a'^ exp J ^ j (35) 
Rescaling a and x in fl35l) yields 

aa(-) + i (5a' + a) a'" = a~' exp J ^ j 
when a 7^ and 

aa^ + -a a = a (37) 

when a = 0. 

When Co 7^ and Ci 7^ 0, we apply transformation ( 1^ to and we obtain 
(again with the rescaling of A and X) 

AA^'-^ + i (5A' - «) A'" = A~^ exp | ^ j (38) 

when a 7^ and 

2 

when a = 0. Changing the sign of a in (!38l) brings it to the form of (!36|) . and so the 
critical equations are and (1571) . 



(36) 



Noether and Lie symmetries for charged perfect fluids 



10 



5.1.1. Case la: a = 0. If a = then (1311) becomes 

Y" + MY = /sF^ + g2Y^ + N (39) 

The solution of fl39l) is expressed as the quadrature 

r dY 
X-Xo= / , (40) 
^ ^y2Y^ + 1/2^3 - MF2 + 2NY + L 

where L is an arbitrary constant introduced in the first integration of (|39|1 . A full 
discussion of the evaluation of this quadrature can be found in Kweyama et al [21]. If 
we invoke the same transformations in obtaining (!39|) to the first integral (12^ we obtain 

-L = -g^Y^ + -f^Y'' - -MY^ + NY - -Y'^ 
4^ 3'^ 2 2 

which is the intermediate integration between ( l39l) and ( 140|) . This again indicates that 

the Noether results are a subset of the Lie results. 

We now have to determine a{x). We make the observation that if a = and 

go = 165(2 and gi = g2/Cf, (IE]) and ([22]) become ([35]) which reduces to ([37]). We 

therefore consider (!37l) for further analysis and reduction to quadrature. In carrying out 

the Lie analysis of (l37)l . using PROGRAM LIE [27], we find that it has two Lie point 

symmetries, namely 

dx 

d A d 
^2 = X— + -a— 
ox 5 oa 

Usually, when an nth order equation admits an m < n dimensional Lie algebra of 
symmetries, there is little hope for the solution of the equation via those symmetries. 
However, in this case we are able to reduce the equation due to the presence of hidden 
symmetries |23j . 

The symmetry Gi determines the variables for reduction 

u = a, V = a' 

and the reduced equation is 

uW + AuWv" + -u'v^v" + u\v'^ + -uW^ -1 = 

2 2 

This equation admits the following two symmetries 

U, = u—- -V— 

du 4 dv 

2d d 

ou ov 

The variables for reduction via U2 are 

_i 3 , 1 1 

r = u 2t> s = u^v u^v 

2 

and the reduced equation is 

r^s'^s" + r^ss'^ + Ar^s'^s' + rs^ - 1 = (41) 
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The Lie symmetry analysis of f l4ip yields the following two symmetries 
^ d s d 

^ dr 3ds 
^ d s d 

The reduction variables generated by X2 are 

p = rs q = rs' + s 
and the reduced equation is 



p^qq' + pq'^ — 1 



with solution 

2 2 go 
p p'^ 

where go is a constant. We can now invert these transformations to find the solution of 
(l37|l . Alternatively, we can integrate (!37|) directly and write down the solution as 

Ux =a-^'^ = [G\u)]-^ 
X — xq = G{u) 

where we have set 

^^^^ ~ / (/C4 + /C3M + /C2M2 - (1/6)/Cim3 + {l/Z2)lCou^f'^ 
and the /Cj, z = 0, . . . , 3 are constants of integration related to M, A^, /2 and g2 and /C4 
is arbitrary. As pointed out earlier, this result was obtained in Kweyama et al [21], but 
for constant c{x). In the case of nonconstant c(x), the solution is the same, except that 
we replace a and x in this solution with A and X respectively. To obtain the solution 
to fl29|) (with a = 0) we need to apply the inverse of fl2T|) . 

5.1.2. Case lb: a ^ 0. When a 7^ 0, we cannot directly reduce fl3Tl) to quadratures. 
We need to investigate the constraints under which it possesses a second point symmetry. 
We find that if /2 7^ 0, 7^ then (13T|) has the following two symmetries 

G. - ^ (42«) 

G, = eW3)^ A _ eW3)^ (^y + ^]^ (426) 
dX \3 9g2j dY ^ ^ 

provided the following conditions are satisfied 

M=-Il-^, N=JL-^J^ (43) 
3g2 36' 21 gl 27g2 ^ ' 

Utilising (142 h\j we obtain the transformation 

X = - ie-("/3)^, y = e^"/^')^ {y + ^\ 
a \ ^92 
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and equation (pTj) becomes 

y" = g2y' 

with solution 

dy 







The values in fH3|) look rather peculiar. However, it is interesting to note that these 
values correspond directly to a simplification of the eigenvalue problem associated with 
a dynamical systems analysis of (1311) ! 

When c 7^ 0, = 0, /2 7^ we find that g = Q (and so we are in the neutral perfect 
fiuid realm). Now (131 p has the following two symmetries 

(44a) 

' dX^ 1^500/2 5/2 5 J dV ^ ^ 

provided the following condition is satisfied 

N = — + —— + 



4/2 8/2 40000/2 
36a^ 




(45) 

2500/2 

This condition is equivalent to the one obtained by Mellin et al ^24j for the case where 
n = 2 in their analysis of the generalised Emden-Fowler equation. We use fl44 h\\ to 
obtain the following transformation 

' ^ I 2/2 2oo/2r 

which, together with ( H5l) reduces equation ( l3T|) to 

y = /23^' 



with solution 



dy 



^y^ + Ci 



We take this opportunity to make two minor corrections to the work of Mellin et al 
|24j : While the expression for given by their equation (7.11) is correct, it is obtained 
by multiplying their equation (7.9) by a / (ia~^/^exp \ f {p — 2Co/a) dx dx and then 
integrating, not multiplying by a / da~^^'^dx as indicated in their paper. Also, the 
coefficient of Co should be 2 in their equation (7.8). 
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5.2. Case II: c = 



From (I25ajl . (!25l . (!25dll and ( l25ell we have 
1 
2 



6 = i (a' + tt) (46a) 



a = oq + ctix + a2a; (466) 
^ =^2a~=^exp (- ) (46cO 



a J 

The symmetry f|T8l) now takes the form 

„ 9 1 , , ,9 
G = a— + - (a + «)?/— 
ax 2 a|/ 

Using the transformation 

adx^ 
~2a 

equation ([6]) is transformed into the autonomous equation 

Y" + aY' + /3Y = f^Y^ + g^Y^ (47) 

where 

/3 = - (a^ - a) , A = aj- 4aoa2 

In carrying out the standard Lie point symmetry analysis on ( H71) we have the following 
cases: 

5.2.1. Case Ila. If /2 7^ 0, (72 7^ 0, fH7|) has the following two symmetries 

G. = A (48a) 

provided the following conditions apply 

4»2 f2 

We use ( 148 fel) to obtain the following transformation 

X = - ^e-W3)X^ y ^ ^-ia/3)X ( Y + ^] (50) 



Using (H9!) and (150|) the equation fj471) becomes 
2a 

and the solution of (15T]) is 



y = ^3^^ (51) 
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5.2.2. Case lib. If f^ = 0, g2 ^ (which imphes that / = 0), then (HZD has the 
following two symmetries 

= ^ (520) 

G. = e<"/3)-^A_2:e<°">-^y 9 (52(,) 
oX 3 ar 

subject to the following condition 

P = (53) 



Using fl526p we obtain the following transformation 



;f = --e-("/3)^, 3; = e("/3)^F (54) 
a 

Using (l53i) and (15^ equation f j471) is transformed to 

y = 92y^ (55) 

and the solution of (155|) is 

X — Xq = ( — , = 

This is an intrinsically charged result - there is no uncharged analogue. 

5.2.3. Case lie. If /2 7^ 0, 5^2 = (which implies 5^ = 0), then (H71) has two symmetries 
provided 

,3 = ±^ (56) 
and can be transformed to 

y" = f2y' 

which has the solution 



X — X( 







dy 



^y^ + c 



This result was previously obtained in the neutral case by Maharaj et al [12]. A 
consequence of (156|1 is that A > (Note that this is not imposed on (147|) as was done in 
Maharaj et al [IS]-) and hence a in (j466j) has real roots. 



6. Discussion 



We have undertaken a comprehensive analysis of in order to determine which forms 
of the functions / and g would lead to first integrals and/or solutions of the equation. We 
reviewed our previous ad hoc approach, and showed that those results were contained in 
the results obtained via Noether's theorem. This occurred when the function c obtained 
in the Noether analysis was set to j. These latter results were then shown to be further 
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contained with those obtained via the Lie analysis. This occurs when a = 0. In this 
case (EZD and ([28]) take the forms of ffT96j) and ffT9c|l respectively, and so (JUft]) and f lT9c|l 
are special cases of ( 127|) and f l28l) respectively. 

Usefully, the first two approaches, namely the ad hoc approach and the Noether 
symmetry analysis, yielded first integrals directly. The final (Lie) approach needed two 
stages of analysis in order to reduce the equation to quadratures. 

We were able to completely analyse ([6]) for Lie point symmetries in an exhaustive 
analysis. All possible cases for the functions / and g were analysed. We showed that 
these results reduced to the neutral case results {g = 0) of Maharaj et al [19] in most 
cases. However, we were also able to find an inherently charged case in §5.2.21 that has 
no uncharged analogue. 

While a complete analysis was produced for the case of a = 0, only a partial 
analysis could be performed when a 7^ 0. Nonetheless, we were still able to provide 
conditions under which ^ could be reduced to quadratures. It still remains to solve 
( l36l) . The main difficulty is that this equation in an integro-differential equation and 
such equations are notoriously difficult to solve. However, we can still reduce fl36]) to a 
first order equation by eliminating a" from ( l32l) and (133]) . Further work in this direction 
is ongoing. 

Note that, in the case that c = 0, we were able to find constraints under which we 
could reduce (16]) to quadratures with / and g given explicitly and a 7^ 0. It is interesting 
to observe that the constraints we found in all subcases of §5.21 forced the quadratic a 
in ( 14661) to have real roots. 

To complete our analysis, we take the symmetries calculated in §5.1.11 and 
investigate the possibility of group invariant solutions [25] [26] of ( l37j) . The only two 
results of significance are that 

^-^.-^/ , , ,,3 w 

_f|y a-(i/3) + ir,a2/3 



and 



x^/^ (58) 



are invariant solutions of (1371) . The first solution can only be given implicitly and so is 
not of much practical use. However, if we let c = | in ([26]) . a = in ( [27]) and Ki = Q 
in 057]), then we obtain explicit forms of / and g. Also if we substitute ([S5[) into ([2B[) 
and ( [27[) we can explicitly determine / and g. In both cases we obtain the forms 

f{x) cx (a; — xqY^^^^^'' , g{x) oc (a; — xq)"*'^^'^^'' (59) 

These forms for / and g were earlier obtained by Kweyama et al [2T] using an ad 
hoc approach which yielded a new charged first integral to the Einstein-Maxwell field 
equations. In this case, ([6]) admits two Lie point symmetries and can be reduced to 
quadratures. 
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We can also reuse (13^ by invoking (^T]) to obtain 

-(6/5) 

(60) 




which is a solution to (129|) provided 

24/?^ , , 

= -625CF- ^^^^ 
This corresponds to the choice 

f{x) (X {x — xq)~^^^^^^ , g{x) oc (x — Xo)"*'"'^'^^^'' (62) 

for which again has two Lie point symmetries and can be reduced to quadratures. 

In summary, we have given a complete Noether and Lie point symmetry analysis 
of (j6]). For the Lagrangian ( [T71) we were able to provide the most general Noether 
point symmetry, the most general first integral associated with this symmetry and 
indicated that this integral was equivalent to that found via an ad hoc approach, i.e 
f ll3p . Finally we determined the most general Lie point symmetry admitted by 
and gave conditions under which the equation could be reduced either to a first order 
equation, or to quadratures. 
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